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Q-f Using dynamic paraboloid model of Earth's magnetosphere, a large set of particles' trajectory 

■ computations was performed. Based on its result, the numerical algorithm for calculating 
effective cutoff rigidity dependence on empiric parameters has been developed for further use in 

Q H " magnetospheric transmission calculations. 
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lB? 1 Introduction 



^ . Radiation condition onboard LEO spacecrafts is determined, particularly, by the charged particles 
CN penetrating from outside of the magnetosphere. The measure of such a process in any given near- 
^ Earth space location is a local geomagnetic cutoff rigidity value. In this work we assume, that 
CO ' "cutoff rigidity" phrase denotes an effective vertical cutoff rigidity (i2 e //), which value is gained from 
VO calculated discrete penumbra structure by ordinary technique using white spectra (for example, [1]). 

It is well known, that magnitude of R e ff an d penetration boundaries' position depend on local 
O , time and different magnetospheric condition-related parameters [21 El IH El El El CD] • Cutoff rigidity 
value's variations can be measured in experiments or obtained by some type of computation using 
Earth's magnetosphere model. However, usual method for R e ff calculation is based on resource 
consuming trajectory computation technique [3 El- 
ba this work we offer the method to calculate R e ff with accounting for local time and main 
empirical parameters, that characterizes magnetospheric condition, in any point of near-Earth's 
space. For the applications where R e ff calculation is often needed (e.g. transmissions for LEOs) 
our method provides a huge speedup in comparison with direct trajectory computations. The only 
things needed for the method are the IGRF rigidity in exploring point and rigidity attenuation [12] 
quotient's dependence on empirical parameters. The first one can be obtained using interpolation 
in Rif^ tabulated data (see, for example, (THUS]), which is regularly updated for coming epochs. 
And the second can be obtained only once by using particle trajectory computations in corresponding 
numerical magnetospheric model. In our case the dynamic paraboloid model [161 El HE] of Earth's 
magnetosphere was chosen for solution of the problem. 
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2 Cutoff rigidity calculation scheme 



We have adopted in our work the cutoff rigidity attenuation formalism [12j [13] , which has already 
successfully applied [15] for transmission calculations using Tsyganenko-89 [14J model. The rigidity 
attenuation quotient A in our case is calculated as 

A modei = Kffifaia)/Bffi (disturbed) 

where R^ odel (quiet) value corresponds the set of model parameters for non-disturbed magnetosphere 
condition, and R^ odel (disturbed) is correspond to some set of varied parameters regarding to quiet 
ones. 

As it was wordlessly postulated in [12], the A dependence on point's geographical position might 
not be taken into account because it is relatively small. Generally, it seems to be true at least for 
R > 0.6 GV with good enough accuracy (error's order is about percents even for lowest R in this 
range). Hence, the approach of calculating a world-wide R e ff grid can be essentially simplified. 

The algorithm we propose for obtaining R^ odel (LT, parameters) in given point contains 3 stages: 

• Calculating Riff 1 ^ for the given point; for practical needs, the interpolation in pre-calculated 
grid for some altitude Ho is often used, it is especially easy when large amount of results is 
needed quickly. The formula 

tdIGRF i TT\ _ tjIGRF i tt \ ( r E + 

R eff ( H ) - Kff i^oj ■ y rE + H J 

can be applied to transform R^f F value to new altitude H. 

• Calculating a model-dependent cutoff rigidity attenuation quotient 
A WO(M (i^gF,LT, parameters) = R™f f de ' '(quiet) / 'R™ff ie ' l (R[ c f '/ lF \LT, parameters) by interpola- 
tion in pre-calculated R™ dd database (or by extrapolation if Rifp F value is out of applied 
basic R^ff F array, see Appendix) 

• Computing the resulting value R™j el = R I Jfj lF / A model (R 1 ^^ , LT, parameter s) 

Table of Rifp F values for some given altitude Hq is needed to be renewed every 5 years with new 
IGRF epoch coming. Values R^j odel (quiet) were obtained by trajectory computations using dynamic 
paraboloid model with parameter set [7, 380, 0, 0, 5, —3, 1] (sequence of variables is according to tabled] 
below). To compute R^j el database which is intended to provide dependence of A model on parameters 
and local time LT, we used the varying of all model parameters in wide enough ranges, comprising 
extreme, quiet and "anti-extreme" parameter sets. The multidimensional grid, applied for this 
calculation, is presented in table [H Values at the edges of this grid were obtained (approximately) 



Table 1: The applied parameter grid. 
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from [19]. All calculations were performed for local time values 3,9,15 and 21 hours. Technical 
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Figure 1: Partial dependencies of R e ff on model parameters for six explored points. 
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details of calculations are summarized in Appendix. Let us now demonstrate some results of cutoff 
rigidity dependence on local time and model parameters, that have been computed using presented 
technique. 

Fig{T] gives an opportunity to weigh deposit of every model parameter in resulting R e ff value for 
all points of noted above set. Such case of partial dependencies is obtained by fixating all parameters 
as quiet and excepting the given one. Note, that the role of Dst as R e ff depression factor wanes with 
decreasing of R e ff, but the AL one rises rapidly. The weight of B IMF and B IMF seems to be not 
important for R e ff values, at any hand, in presented case of quiet basic values, because they result 
in R e ff changes that are comparable with the value of applied rigidity step size (see Appendix). 
Contrary, decrease of B IMF effects in huge R e ff grow (see also fig 12). The effect of p and v here is 
clear but of moderate magnitude. 

Fig f2] shows family of R e ff daily variations for three points having normal cutoff values of different 
order. For the first picture there are some hours where R e ff is negligibly small, so such a point can 
be accessible for particles of any rigidity. The middle picture demonstrates the effect of B IMF 
"anti-extreme" setting, leading not only to very high R e ff mean values (for this case), but also to 
the translocation of curve's maximum for time axis's positive direction. The upper curve on the 
most right part of FigfS] exhibits the insufficiency of applied step size for penumbra calculation (see 
Appendix), leading to the underestimation of R e ff value for LT = 9, which causes an ambiguous 
interpolation results. However, this effect is rare and too weak to affect resulting transmission, for 
which calculating presented method was created. 




Local Time, hours Local Time, hours Local Time, hours 

Figure 2: Dependencies of R e ff on local time and model parameter set for three points with 
R™f f del (quiet) values 0.49 GV (left), 1.98 GV (middle) and 5.48 GV (right). 

A "gnarlyness" of some presented curves is due to effect of scarce grid, currently imperfect 
interpolation procedure and scale. It is does not sufficiently affect the results of applications where 
much of R e ff values to be calculated, because interpolation uses many R^ odel values and their errors 
(which are of order of R step size) are to be evened. Extrapolating parameter values to out of grid 
(but not very far) is also possible and seems to be correct enough to get reliable result. 
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3 Transmission functions for LEO missions 



Magnetospheric transmission functions are directly connected with such tasks as LEO radiation 
conditions evaluating and the interpretation of orbital experiment results [21]. Because of their 
importance, it is necessary to be able to account for more empiric factors that affecting transmission. 
Here we'll give some examples of calculated magnetospheric transmission functions for several (quiet, 
disturbed and very extreme) model parameters sets for ISS-like orbit, obtained using presented 
technique. The essential note is that, unlike Tsyganenko models, the paraboloid one is primordially 
intended to describe even very extreme conditions, as it is not internally limited by experimental data 
arrays. Figj3] presents some modeled transmissions for H = 450 km i = 51.7° circular orbit under 
different conditions for long mission duration, where all localities are averaged and transmission curve 
became smoothed. In comparison with additionally figured transmission, obtained by Tsyganenko-89 
(T89) model with Kp = 6, it is obviously seen that the paraboloid model gives unapproachable for 
T89 results in modeling most extreme conditions. 




Figure 3: Magnetospheric transmission ^(R) for different conditions, obtained according to 
paraboloid model with using presented technique. Tsyganenko-89 transmission for Kp = 6 is also 
figured here for comparison. Note also the standing out of all " anti-extreme" case. 



4 Conclusion 

The method we have developing is intended for applications using intensive geomagnetic cutoff calcu- 
lations, first of all for express magnetospheric transmission calculation for given LEO missions, even 
under changing magnetospheric conditions during flight. Test R e ff and transmission calculation, 
based on the presented technique, is available online by URL |http : / / dec 1 . sinp . msu . ru/ ~ vovka / riho 
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(a simplified version, where transmission calculation is allowed for only static conditions and circular 
orbits). 

Current version of the method uses the IGRF rigidity values table as basic (in step 1 of scheme), 
although it seems that the best choice would be to apply some other world basic grid, for example, 
one directly calculated R™^ 1 (quiet) table with using paraboloid model. Nevertheless, it does not 
understate the attenuation quotient methodology itself, because its drawback in our case is only 
correspond to relatively small systematic inaccuracies between conditional " quiet" R™j el values and 
IGRF ones. 
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Appendix. Technical details 

Here we summarize technical moments that are related to the performed R^°j^ el trajectory calculations 
with using paraboloid model. 

Vertically directed protons were test particles for reverse trajectory calculations of penumbra. 
Fourth order integration scheme was used for it, with accounting for time changing (hence, magnetic 
field) during particle's modeled flight. The geomagnetic field there was superposition of IGRF epoch 

2005 field with dynamic paraboloid model (version 2004 [20]), with uniform field vector B IMp outside 
of the magnetopause, which position is natively given by the code realizing the paraboloid model. 
Maximal flying distance was equal to 15rg, the particles walked it over during motion were considered 
as allowed for penetration. The Earth was represented by WGS-82 ellipsoid with atmosphere layer 
at 20 km above its surface, the particles which fell below it was considered as reentrant. Rigidity 
step size for penumbra calculation was equal to 0.01 GV. All calculations were performed for altitude 
450 km above mean Earth's radii = 6371.2 km for six selected points in northern geographic 
hemisphere with basic IGRF rigidities 12.086,7.381,4.101,2.018,1.203,0.666 GV. 

All presented transmissions were obtained for mission duration 3000 revolutions with orbital 
trajectory step size 0.9 degree. 
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